Solar Multi-Scale Convection and Rotation Gradients Studied in Shallow Spherical Shells

The differential rotation of the sun, as deduced from helioseismology, exhibits a prominent radial shear layer near the top of the convection zone wherein negative radial gradients of angular velocity are evident in the low- and mid-latitude regions spanning the outer 5% of the solar radius. Supergranulation and related scales of turbulent convection are likely to play a significant role in the maintenance of such radial gradients, and may influence dynamics on a global scale in ways that are not yet understood. To investigate such dynamics, we have constructed a series of three-dimensional numerical simulations of turbulent compressible convection within spherical shells, dealing with shallow domains to make such modeling computationally tractable. These simulations are the first models of solar convection in a spherical geometry that can explicitly resolve both the largest dynamical scales of the system (of order the solar radius) as well as smaller-scale convective overturning motions comparable in size to solar supergranulation (20--40 Mm). We find that convection within these simulations spans a large range of horizontal scales, and that the radial angular velocity gradient in these models is typically negative, especially in low- and mid-latitude regions. Analyses of the angular momentum transport indicates that such gradients are maintained by Reynolds stresses associated with the convection, transporting angular momentum inward to balance the outward transport achieved by viscous diffusion and large-scale flows in the meridional plane. We suggest that similar mechanisms associated with smaller-scale convection in the sun may contribute to the maintenance of the observed radial shear layer located immediately below the solar photosphere.Comment: 45 pages, 17 figures, ApJ in press. A preprint of paper with hi-res figures can be found at http://www-lcd.colorado.edu/~derosa/modelling/modelling.htm

Structure and Evolution of Giant Cells in Global Models of Solar Convection

The global scales of solar convection are studied through three-dimensional simulations of compressible convection carried out in spherical shells of rotating fluid which extend from the base of the convection zone to within 15 Mm of the photosphere. Such modelling at the highest spatial resolution to date allows study of distinctly turbulent convection, revealing that coherent downflow structures associated with giant cells continue to play a significant role in maintaining the strong differential rotation that is achieved. These giant cells at lower latitudes exhibit prograde propagation relative to the mean zonal flow, or differential rotation, that they establish, and retrograde propagation of more isotropic structures with vortical character at mid and high latitudes. The interstices of the downflow networks often possess strong and compact cyclonic flows. The evolving giant-cell downflow systems can be partly masked by the intense smaller scales of convection driven closer to the surface, yet they are likely to be detectable with the helioseismic probing that is now becoming available. Indeed, the meandering streams and varying cellular subsurface flows revealed by helioseismology must be sampling contributions from the giant cells, yet it is difficult to separate out these signals from those attributed to the faster horizontal flows of supergranulation. To aid in such detection, we use our simulations to describe how the properties of giant cells may be expected to vary with depth, how their patterns evolve in time, and analyze the statistical features of correlations within these complex flow fields.Comment: 22 pages, 16 figures (color figures are low res), uses emulateapj.cls Latex class file, Results shown during a Press release at the AAS meeting in June 2007. Submitted to Ap

Inferring Maps of the Sun's Far-side Unsigned Magnetic Flux from Far-side Helioseismic Images using Machine Learning Techniques

Accurate modeling of the Sun's coronal magnetic field and solar wind structures require inputs of the solar global magnetic field, including both the near and far sides, but the Sun's far-side magnetic field cannot be directly observed. However, the Sun's far-side active regions are routinely monitored by helioseismic imaging methods, which only require continuous near-side observations. It is therefore both feasible and useful to estimate the far-side magnetic-flux maps using the far-side helioseismic images despite their relatively low spatial resolution and large uncertainties. In this work, we train two machine-learning models to achieve this goal. The first machine-learning training pairs simultaneous SDO/HMI-observed magnetic-flux maps and SDO/AIA-observed EUV 304$\r{A}$ images, and the resulting model can convert 304$\r{A}$ images into magnetic-flux maps. This model is then applied on the STEREO/EUVI-observed far-side 304$\r{A}$ images, available for about 4.3 years, for the far-side magnetic-flux maps. These EUV-converted magnetic-flux maps are then paired with simultaneous far-side helioseismic images for a second machine-learning training, and the resulting model can convert far-side helioseismic images into magnetic-flux maps. These helioseismically derived far-side magnetic-flux maps, despite their limitations in spatial resolution and accuracy, can be routinely available on a daily basis, providing useful magnetic information on the Sun's far side using only the near-side observations.Comment: Accepted by Ap

A Method for Data-Driven Simulations of Evolving Solar Active Regions

We present a method for performing data-driven simulations of solar active region formation and evolution. The approach is based on magnetofriction, which evolves the induction equation assuming the plasma velocity is proportional to the Lorentz force. The simulations of active region coronal field are driven by temporal sequences of photospheric magnetograms from the Helioseismic Magnetic Imager (HMI) instrument onboard the Solar Dynamics Observatory (SDO). Under certain conditions, the data-driven simulations produce flux ropes that are ejected from the modeled active region due to loss of equilibrium. Following the ejection of flux ropes, we find an enhancement of the photospheric horizontal field near the polarity inversion line. We also present a method for the synthesis of mock coronal images based on a proxy emissivity calculated from the current density distribution in the model. This method yields mock coronal images that are somewhat reminiscent of images of active regions taken by instruments such as SDO's Atmospheric Imaging Assembly (AIA) at extreme ultraviolet wavelengths.Comment: Accepted to ApJ; comments/questions related to this article are welcome via e-mail, even after publicatio

Implications of Different Solar Photospheric Flux-Transport Models for Global Coronal and Heliospheric Modeling

The concept of surface-flux transport (SFT) is commonly used in evolving models of the large-scale solar surface magnetic field. These photospheric models are used to determine the large-scale structure of the overlying coronal magnetic field, as well as to make predictions about the fields and flows that structure the solar wind. We compare predictions from two SFT models for the solar wind, open magnetic field footpoints, and the presence of coronal magnetic null points throughout various phases of a solar activity cycle, focusing on the months of April in even-numbered years between 2012 and 2020, inclusive. We find that there is a solar cycle dependence to each of the metrics considered, but there is not a single phase of the cycle in which all the metrics indicate good agreement between the models. The metrics also reveal large, transient differences between the models when a new active region is rotating into the assimilation window. The evolution of the surface flux is governed by a combination of large scale flows and comparatively small scale motions associated with convection. Because the latter flows evolve rapidly, there are intervals during which their impact on the surface flux can only be characterized in a statistical sense, thus their impact is modeled by introducing a random evolution that reproduces the typical surface flux evolution. We find that the differences between the predicted properties are dominated by differences in the model assumptions and implementation, rather than selection of a particular realization of the random evolution.Comment: Accepted for publication in The Astrophysical Journa

Global MHD Simulations of the Time-Dependent Corona

We describe, test, and apply a technique to incorporate full-sun, surface flux evolution into an MHD model of the global solar corona. Requiring only maps of the evolving surface flux, our method is similar to that of Lionello et al. (2013), but we introduce two ways to correct the electric field at the lower boundary to mitigate spurious currents. We verify the accuracy of our procedures by comparing to a reference simulation, driven with known flows and electric fields. We then present a thermodynamic MHD calculation lasting one solar rotation driven by maps from the magnetic flux evolution model of Schrijver & DeRosa (2003). The dynamic, time-dependent nature of the model corona is illustrated by examining the evolution of the open flux boundaries and forward modeled EUV emission, which evolve in response to surface flows and the emergence and cancellation flux. Although our main goal is to present the method, we briefly investigate the relevance of this evolution to properties of the slow solar wind, examining the mapping of dipped field lines to the topological signatures of the "S-Web" and comparing charge state ratios computed in the time-dependently driven run to a steady state equivalent. Interestingly, we find that driving on its own does not significantly improve the charge states ratios, at least in this modest resolution run that injects minimal helicity. Still, many aspects of the time-dependently driven model cannot be captured with traditional steady-state methods, and such a technique may be particularly relevant for the next generation of solar wind and CME models

Inferring the structure of the solar corona and inner heliosphere during the Maunder minimum using global thermodynamic magnetohydrodynamic simulations

Observations of the Sun’s corona during the space era have led to a picture of relatively constant, but cyclically varying solar output and structure. Longer-term, more indirect measurements, such as from 10Be, coupled by other albeit less reliable contemporaneous reports, however, suggest periods of significant departure from this standard. The Maunder Minimum was one such epoch where: (1) sunspots effectively disappeared for long intervals during a 70 yr period; (2) eclipse observations suggested the distinct lack of a visible K-corona but possible appearance of the F-corona; (3) reports of aurora were notably reduced; and (4) cosmic ray intensities at Earth were inferred to be substantially higher. Using a global thermodynamic MHD model, we have constructed a range of possible coronal configurations for the Maunder Minimum period and compared their predictions with these limited observational constraints. We conclude that the most likely state of the corona during—at least—the later portion of the Maunder Minimum was not merely that of the 2008/2009 solar minimum, as has been suggested recently, but rather a state devoid of any large-scale structure, driven by a photospheric field composed of only ephemeral regions, and likely substantially reduced in strength. Moreover, we suggest that the Sun evolved from a 2008/2009-like configuration at the start of the Maunder Minimum toward an ephemeral-only configuration by the end of it, supporting a prediction that we may be on the cusp of a new grand solar minimum